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Abstract 

The exact solution of the asymmetric exclusion problem with first- and 
scond-class particles is presented. In this model the particles (size 1) of both 
classes are located at lattice points, and diffuse with equal asymmetric rates, 
but particles in the first class do not distinguish those in the second class from 
holes (empty sites). We generalize and solve exactly this model by considering 
molecules in the first and second class with sizes si and ,32 (si, S2 = 0, 1, 2, . . .), 
in units of lattice spacing, respectively. The solution is derived by a Bethe 
ansatz of nested type. We give a simple pedagogical presentation of the Bethe 
ansatz solution of the problem which can easily be followed by a reader not 
specialized in exactly integrable models. 

I. INTRODUCTION 

The similarity between the master equation describing time fluctuations in nonequilib- 
rium problems and the Schrodinger equation describing the quantum fluctuations of quan- 



turn spin chains turns out to be fruitful for both areas of research [jT| - Since many 

quantum chains are known to be exactly integrable through the Bethe ansatz, this provides 
exact information on the related stochastic model. At the same time the classical physical 
intuition and probabilistic methods successfully applied to nonequilibrium systems give new 
insights into the understanding of the physical and algebraic properties of quantum chains. 

An example of this fruitful interchange is the problem of asymmetric diffusion of hard- 
core particles on the one dimensional lattice. This model is related to the exactly integrable 
anisotropic Heisenberg chain ]16|] (XXZ model). However if we demand this quantum chain 
to be invariant under a quantum group symmetry Uq{SU{2)), we have to introduce, for the 
equilibrium statistical system, unphysical surface terms, which on the other hand have a 
nice and simple interpretation for the related nonequilibrium stochastic system 0,^. 

In the area of exactly integrable models it is well known that one of the possible extensions 
of the spin-| XXZ chain to higher spins, is the anisotropic spin-1 Sutherland model (grading 
^1 = ^2 = ^3 = 1) other hand in the area of diffusion limited reactions 

a simple extension of the asymmetric diffusion problem is the problem of diffusion with 



first and second class of particles []TB[ - pO] . In this problem a mixture of two classes of 



hard-core particles diffuses on the lattice. Particles belonging to the first class ignore the 
presence of those in the second class, i. e., they see them in the same way as they see the 
holes (empty sites). In it was shown that for open boundary conditions the anisotropic 
spin-1 Sutherland model and this last stochastic model are exactly related, the Hamiltonian 
governing the quantum or time fluctuations of both models being given in terms of generators 
of a Hecke algebra, invariant under the quantum group UqSU{3). In this paper we are going 
to derive the exact solution of the associated quantum chain, in a closed lattice. Recently 



p!5| (see also [|n|) we have shown that without losing its exact integrability, we can consider 
the problem of asymmetric diffusion with an arbitrary mixture of molecules with different 
sizes (even zero), as long they do not interchange positions, that is, there is no reactions. 
Motivated by these results we are going to extend the asymmetric diffusion problem with a 
second class of particles, to the case where the particles in each class have an arbitrary size. 



in units of the lattice spacing. Unlike the case of asymmetric diffusion problem, we have in 



this nested Bethe ansatz ||2T|. 



The paper is organized as follows. In the next section we introduce the generalized 
asymmetric model with second-class particles and derive the associated quantum chain. In 
section 3 the Bethe ansatz solution of the model is presented in a pedagogical and self 
contained way, which should be simple to follow by an audience not specialized in exactly 
integrable models. Finally in section 4 we present our conclusions, with some possible 
generalizations of the stochastic problem considered in this paper, and some perspectives on 
future work. 

II. THE GENERALIZED ASYMMETRIC DIFFUSION MODEL WITH FIRST- 

AND SECOND-CLASS OF PARTICLES 

A simple extension of the asymmetric exclusion model, where hard-core particles diffuse 
on the lattice, is the problem where a mixture of particles belonging to different classes (first 
and second class) diffuses on the lattice. This problem was used to describe shocks |0 - ||2^ 
in nonequilibrium and also has a stationary probability distribution that can be expressed 
via the matrix-product ansatz . In this model we have rii and n2 molecules belonging to 
the first and second class, respectively. Both classes of molecules diffuse asymmetrically, but 
with the same asymmetrical rates, whenever they encounter empty sites (holes) at nearest- 
neighbor sites. However, when molecules of different classes are at their minimum separation, 
the molecules of the first class exchange position with the same rate as they diffuse, and 
consequently the first-class molecules see no difference between molecules belonging to the 
second class and holes. 

We now introduce a generalization of the above model, where instead of having unit 
size, the molecules in the first and second class have in general distinct sizes si and S2 
(si, S2 = 1,2,...), respectively, in units of lattice spacing. In Fig. 1 we show some examples 
of molecules of different sizes. We may think of a molecule of size s as formed by s monomers 
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(size 1), and for simplicity, we define the position of the molecule as the center of its leftmost 
monomer. The molecules have a hard-core repulsion: the minimum distance dap, in units 
of the lattice spacing, between molecules a and (3, with a in the left, is given by da/s = Sa- 
In order to describe the occupancy of a given configuration of molecules we define, at each 
site z of a lattice with sites, a variable (3i {i = 1,2, ... , N), taking the values Pi = 0, 1 
and 2, representing site i empty (size sq = 1), occupied by a molecule of class 1 (size si) 
or a molecule of class 2 (size S2), respectively. Then the allowed configurations are given 
by the set {i = 1, . . . ,N), where for each pair (/3j, 7^ with j > i we should have 

The time evolution of the probability distribution P{{(3},t), of a given configuration {(3} 
is given by the master equation 

^^^^^ = -r({/^} - {f3'})pm,t)+rm ^ m)pm,t), (i) 

where r({/?} — > is the transition rate for configuration {(3} to change to {/?'}. In 

the present model we only allow, whenever it is possible, the particles to diffuse to nearest- 
neighbor sites, or to exchange positions. The possible motions are diffusion to the right 

li 0i+si 0i li+i, (rate Tn) 

2i 0,+,, -^0i 2,+i, (rate Vr), (2) 

diffusion to the left 

0i li+i li %i+s^, (rate r^) 

0, 2,+i ^2i 0,+,,, (rate r^), (3) 
and interchange of particles 

li 2j_|_sj — >■ 2j Ii+S2) (rate Tr) 

2i Ii+S2^h 2i+si, (rate r^). (4) 

As we see from (^), particles in the first class interchange positions with those of second 
class with the same rate as they interchange positions with the empty sites (diffusion). We 



should remark however that unless the second class particles have unit size (s2 = 1), the 
net effect of the second class particles in those of the first class is distinct from the effect 
produced by the holes, since as the result of the exchange the first class of particles will 
move by S2 lattice size units, accelerating its diffusion. 

The master equation (0) can be written as a Schrodinger equation in Euclidean time (see 
Ref. for general application for two body processes) 

if we interpret |P >= P{{(3},t) as the associated wave function. If we represent (3i as \P >i 
the vector |/3 >i ®\l3 >2 ® ■ ■ ■ ® 1/3 >Ar will give us the associated Hilbert space. The process 
(I)- (I) gives us the Hamiltonian (see Ref. for general applications) 

j 
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with 

D = Tn + TL, e+ = -^^, e_ = (e+ + 6_ = 1), (7) 

and periodic boundary conditions. The matrices are 3x3 matrices with a single 

nonzero element (i?°'^)j,,- = Sa^iSpj {a,P,i,j G Z). The projector V in (H), projects out 
from the associated Hilbert space the vectors \{f3} > which represent forbidden positions 
of the molecules due to their finite size, which mathematically means that for all i,j with 
Pi, [3j 7^ 0, \i — j\ > S[^-{j > The constant D in (H) fixes the time scale and for simplicity 
we chose D = 1. A particular simplification of (^) occurs when molecules of the first and 
second class have the same size Si = S2 = s. In this case the Hamiltonian can be expressed 
as an anisotropic nearest-neighbor interaction spin-1 SU{3) chain. Moreover in the case 
where their sizes are unity (s = 1) the model may be related to the SU{3) Sutherland model 
with twisted boundary conditions fl^ . 



III. THE BETHE ANSATZ EQUATIONS 



We present in this section the exact solution of the general quantum chain (|^). Since the 
present paper is going to appear in a special issue where the main subject is nonequilibrium 
physics, we are going to present a pedagogical and self-contained derivation of the exact 
solution, that can clearly be followed by a nonexpert in the arena of exactly solvable models 
in statistical mechanics. 

Due to the conservation of particles in the diffusion and interchange processes the total 
number of particles rii and n2 of particles of class 1 and 2 are good quantum numbers and 
consequently we can split the associated Hilbert space into block disjoint sectors labeled by 
the numbers ni and n2 {fii = 0, 1, . . . , n; n2 = n — ni, n = ni + ■ We therefore consider 
the eigenvalue equation 

H\ni,n2 >= E\ni,n2 >, (8) 

where 

\ni,n2 >= ^^/(xi,Qi;...;a:„,(5„)|xi,(5i;...;x„,(5„ > . (9) 

{Q} M 

Here \xi,Qi; . . . ; Xn,Qn > means the configuration where a particle of class Qi{Qi = 1,2) 
is at position Xi [xi = 1,...,N). The summation {Q} = {Qi, . . . ,Qn} extends over all 
permutations of n numbers in which rii terms are 1 and n2 terms are 2, while the summation 
{x} = {xi, . . . , Xn} runs, for each permutation {Q}, in the set of the n nondecreasing integers 
satisfying 

Xi+i>Xi + sq^, i = l,...,n-l, Sq^ < Xn - Xi < N - Sq^. (10) 

Before getting the results for general values of n let us consider initially the cases where we 
have 1 or 2 particles. 

n = 1. For one particle on the chain (class 1 or 2), as a consequence of the translational 
invariance of (^) it is simple to verify directly that the eigenf unctions are the momentum-fc 
eigenf unctions 
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N N 

\l,0>=J2f{x,l)\x,l>, or |0,l>=^/(x,2)|x,2>, (11) 

x=l x=l 

with 

2ttI 

/(x,l) = /(x,2) = e^'=^ fc=— ,/ = 0,l,...,iV-l, (12) 
and energy given by 

E = e{k) = -{e-e''' + e+e-''' -1). (13) 

n =2. For two particles of classes Qi and Q2 {Qi, Q2 = 1,2) on the lattice, the eigenvalue 
equation gives us two distinct relations depending on the relative location of the particles. 
The first relation applies to the case in which a particle of class Qi (size sqJ is at position 
Xi and a particle Q2 (size sq^) is at position X2, where X2 > xi + sqj. We obtain in this case 
the relation 

Efixi, Qi; X2, Q2) = -e+fixi - 1, Qi, X2, Q2) - e_/(xi, Qi; X2 + 1, Q2) 
- e-/(xi + 1, Qi; X2, Q2) - e+/(a^i, Qi\ - 1, Q2) + 2/(xi, Qi; Xa, Q2), (14) 

where we have used the relation e+ + e_ = 1. This last equation can be solved promptly by 
the ansatz 

f{xi,QuX2,Q2)=e'''''''e'^^^\ (15) 

with energy 

E = e{k{) + e{k2). (16) 

Since this relation is symmetric under the interchange of ki and ^2? we can write a more 
general solution of ([I^ ) as 



/(xi, Qi; X2, g2) = E A«^;g^e*('=^i-^+^'^ 
p 

_ y^Qi,Q2 ^i(kiXi+k2X2) _j_ y^Qi,Q2 ^i{k2Xi+kiX2) (17) 



with the same energy as in (|T6]). In ([I7D the summation is over the permutations P = Pi, P2 



of(l,2). The second relation apphes when X2 = Xi + sq^. In this case instead of (|T^ we have 



Ef{xi,Qi ; xi + sqi, Q2) = -e+f{xi - 1, Qi, xi + sq^,Q2) - e_/(xi, Qi, xi + sq^ + 1, Q2) 
- ^Qi,Q2f{xi^Q2;xi + sq2,Qi) + (1 + eQ,^Q,)f{xi,Qi;xi + sq„Q2), (18) 

where 



61,1 = 62,2 = 0, ei,2 = e+ and £2,1 = e_. 



(19) 



If we now substitute the ansatz (|T7D with the energy (0), the constants Af2'^^ and A^^^^'*^^, 
initially arbitrary, should now satisfy 

E{ l^PuP2 + e'"^- (1 - e~Q„Qj] e^^-.(^«i-i)Agvg^ + ~eQ„Q,e''-^'^^ AfJ^} = (20) 



where 



6+ + e_e' 



i(ki+kj)\ 



(21) 



At this point it is convenient to consider separately the case where Qi = Q2 from those 
where ^ g2. If Qi = <52 = Q (Q = 1, 2) eq. (0) gives 



(22) 



and the cases Qi 7^ Q2 give us the equations 



E 



pikp (si-1) /1 1,2 



jfcp (S2-1) 42,1 



0. 



Performing the above summation we obtain, after lengthy but straightforward algebra, the 
following relation among the amplitudes 



^l,2gifc2(,,_l) 
^2,lgifc2(,2-l) 



Pi,2 + e'"^ 
Pi 2 + e^^^ 



l-$(A;i,A;2 



— e_ e_ 



4I2 jA;i(si-l) 
421 pjfci(s2-l) 



where 



Equations (p2]) and (ESf) can be written in a compact form 



m,h) = ^ (23) 



= -Hp„p, ^ Sg;;g,^(A;p,,A:pj4J;?;\ (gi,Q2 = l,2) (24) 

with 



(25) 



where we have introduced the S matrix. From ^ and (^), the only non zero elements of 



qQl,Q2 



(A:i, k2) = e^('=i-^'2)(«i-i), S2;2(fci, A;2) = e'^'''-'''^^''-^\ 
Sl'Jih, k2) = [1 - e+$(A:i, ^2)] e^('=^-^'^)(^^-i), 
'^1:2(^:1, A:2) = [1 - e^^h,h)] e^(fci-'==^)(^2-i), 
'5i':2(^i,^2) = e+$(A;i,A;2)e^'=^(^^-^)e-^'=^(^^-i), 

S^lih, k2) = e_$(A;i, A;2)e''^^(^^-^)e-'^^(^^-^\ (26) 



A graphical representation of the S matrix is shown in Fig. 2a. Equations (24) do not fix 

the "wave numbers" ki and ^2- In general, these numbers are complex, and are fixed due to 
the cyclic boundary condition 

fixi, Qi; X2, Q2) = /(X2, Q2; xi + N, Qi), (27) 
which from (^) gives the relation 

AQ1Q2 _ ^ikiN aQ2,Qi aQi,Q2 _ Ak2N aQ2,Q\ fOQ\ 

^1,2 — S ^2,1 1 ^2,1 — S ^2,1 • {^°) 



This last equation, when solved by exploiting (P1D-(PBD, gives us the possible values of ki 
and k2, and from (|T^ the eigenenergies in the sector with 2 particles. Instead of solving 
these equations for the particular case = 2 let us now consider the case of general n. 
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General n. The above calculation can be generalized for arbitrary values of ni and n2 
of particles of classes 1 and 2, respectively {rii + n2 = n). The ansatz for the wave function 
(y) becomes 

fix,, Q,;...;xn,Qn) = Yl A^i:;^y^'-^^^+-+'--^-\ (29) 

p 

where the sum extends over all permutations P of the integers 1,2, ... ,n. For the compo- 
nents \xi, Qi, . . . ; Xn, Qn > where Xj+i — Xj > Sq^ for i = 1, 2, . . . , n, it is simple to see that 
the eigenvalue equation (||) is satisfied by the ansatz (p9D with energy 

n 

E = J2eik,). (30) 

On the other hand if a pair of particles of class Qi, Qi+i is at positions Xi, Xi+i, where Xj+i = 
Xi + sq., equation (|) with the ansatz (^) and the relation ( PP] ) give us the generalization 
of relation (P^, namely 

2 

A---,Qi,Qi+i,--- _ cQi'Qi+i/i, h. \ /\---^Q2'Q\^" n D — 1 9 ('^^\ 

^...,Pi,P,+ i,... - --P,,P,+i "^Q'l.Q^ I'^^i'z' '^^'z+J^...,P,+ i,P„... Vi, Vi+i - -L,^, loJ-j 

^1 '^2 



with S given by eq. (^). Inserting the ansatz ( pUj ) in the boundary condition 

/(xi, Qi; . . . ; Q„) = /(xs, Q2; • • • ; a^n, Qn\ Xi + A^, Qi) (32) 
we obtain the additional relation 

/iQl, ■■■,Qn ^ikp.N AQ2,---,Qn,Ql ^'QQ^ 

^Pi,...,p„ - e ^P2,...,p„,Pi > i-^-^; 

which together with (pTD should give us the energies. 

Successive applications of (|3TD give us in general distinct relations between the am- 
plitudes. For example A]]'kf kick's relate to ^^Z^lfk^'ki ^y performing the permutations 
a/37 ~^ ~^ ~^ lf^C( or a/37 ~^ '^iP ~^ 70^/3 — > 7/3a, and consequently the S'-matrix 
should satisfy the Yang-Baxter P^J23[ ] equation 

J2 S"'^^' ih, k2)Sj'"„ {ki, k3)S]j,']j„{k2, ks) = 

7,7',7"=1 

E s^:;^::ik2,k,)s^;j:ik,,ks)sj';';{k,,k2), m 

7,7',7"=1 
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for a,a' , a" , P, (3' , P" = 1,2 and S given by (]26|). In Fig. 2b we show graphically this 
equation. Actually the relation ( plD is a necessary and sufficient condition |p!6| , p3| to obtain 
a non-trivial solution for the amplitudes in Eq. (|3T|). 

We can verify by a long and straightforward calculation that for arbitrary values of Si and 
S2, the S matrix (^Bj), satisfies the Yang-Baxter equation (0), and consequently we may use 
relations (|3l|) and (|33| ) to obtain the eigenenergies of the Hamiltonian (^). Applying relation 



(pl|) n times on the right of equation (^) we obtain a relation between the amplitudes with 
the same momenta. Using the graphical representation in Fig. 3a for the S matrix, we 
illustrate in Fig. 4 the result of such applications. We then obtain 

aQi,-,Q„ _ Akp^N AQ2,-,Qn,Qi _ ( TT \ Jkp^N Sr^ \^ 

Vi=2 / q;,...,q;,q'/,...,q^ 

^Q[,Q'li^Pi^^Pi)^Ql',Q'^i^P2y^Pi) ■ ■ ■ SQ"J^^^^l_^{kp^_^,kp^)SQ7^ (35) 

where we have introduced the harmless extra sum 

2 2 „ 

In order to fix the values of {kj} we should solve (|35|) , i.e., we should find the eigenvalues 
A(/c) of the matrix 

nk)W} = E ( fn SQ^^'ikp,,k^ s$:^{kp„, k)] , (3?) 

Q'l,...,Q'^=l I V«=l / ) 

and the Bethe-ansatz equations which fix the set {ki} will be given from (|35D by 

e-^^^-^ = (-1)"-^ {fl^i,?j A{k,), J = l,...,n. (38) 

Using the graphical representation for S given in Fig. 3a the matrix T in { ^7\j can be 
represented graphically as in Fig. 5, and we identify T{k) as the transfer matrix of an 
inhomogeneous 6- vertex model, on a periodic lattice, with Boltzmann weights S^}''^?{kp^^ k) 
{I = 1, . . . ,n). Consequently, in order to obtain the eigenenergies of the quantum chain (^ 
we should diagonalize the above transfer matrix T[k). 
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Diagonalization of T{k) 

The simplest way to diagonalize T is through the introduction of the monodromy matrix 
Ai{k) |21]], which is a transfer matrix of the inhomogeneous vertex model under considera- 
tion, where, instead of being periodic, the first and last link in the horizontal direction are 
fixed to the values /ii and fin+i (Aii,/^n+i = 1,2), that is 

(39) 



The monodromy matrix ■M.^glj^i^^^ik) is represented in Fig. 6a and has coordinates 
{Q},{Q'} in the vertical space (2" dimensions) and coordinates /ii,/i„+i in the horizon- 
tal space (4 dimensions). This matrix satisfies the following important relations 

This relation is shown graphically in Fig. 6b; its validity follows directly from successive 
applications of the Yang-Baxter equations (Q) (see Fig. 2b). 

In order to exploit relation (HO) let us denote the components of the monodromy matrix 



in the horizontal space by 

m = A(fc)|2>j = A^|:'/}fi(fc), B{k) = 5(fc){j> = M\^:}^i{k), 

c{k) = c{k)\^j\ = m = D(fc)|2>j = M\j:};i{k), (41) 

and clearly the transfer matrix T{k) of the periodic inhomogeneous lattice, we want to 
diagonalize, is given by 

T{k) = A{k) + D{k). (42) 

As a consequence of ( |^ ) the matrices A, B, C and D in ( ^Tf ) obey some algebraic relations. 
The elements (i^i, /ii, Vn+i, fJ^n+i) = (1, 1, 2, 2), (1, 1, 1, 2) and (1, 2, 2, 2) give us the relations 

[S(^),S(A;')]=0, (43) 
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A{k)B{¥) = ||Mls(A:')A(^) - ^^Bik)Aik'), (44) 

^^^)B{k') = ^^Bik')Dik) - ^$^Bik)Dik% (45) 

respectively. The diagonalization of T{k) in ( ^21) will be done by exploiting the above 
relations. This procedure is known in the literature as the algebraic Bethe ansatz |^ . The 



first step in this method follows from the identification of a reference state >, which should 
be an eigenstate of A{k) and D{k), and hence T{k), but not of B{k). In our problem it is 
simple to see that a possible reference state is \VL >= \{ai = 1} ,,,n, which corresponds 
to a state with first-class particles only. It is simple to calculate 

A{k)\n> =a{k)\n>, D{k)\n>= d{k)\n>, 

n 

Cik)\n> =0, B{k)\n>=J2bi{k)\^>. (46) 



i=l 



where 

n n 

<k) = n sij{kp^,k), d{k) = n si:i{kp^,k), 

i=l i=l 
i—1 n 

Hk) = n Slj{kp,,k) n Sl:^ikp^,k). (47) 
1=1 l=j 

In Fig. 7 we show a graphical representation of the above relation. The idea in the algebraic 
Bethe ansatz is that B{k) acts as a creation operator in the reference ("vacuum") state, i.e., 
it creates particles of second class in a sea of particles of first class >, and it is expected 
that a general eigenvector of T{k) in the sector with n2 second-class particles will be given 
by the ansatz 

vl/(fc«, 4^), . . . , kS) = Bik['^)Biki'^) . . . Bik^m > ■ (48) 

The numbers {kl:^^},i = l,...,n2 here play the role of the "wave number" {ki} in the 
coordinate Bethe ansatz (|29|), and are going to be fixed by the eigenvalue equation 

r(fc)vl/(A:«, . . . , kS) = {A{k) + D{k)) ^>{k'l\ . . . , fc«) = A(fc)vl/(fc«, . . . , A;«). (49) 
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Before deriving the relation for general values of n2, let us consider initially the cases 
71-2 = 1 and n2 = 2. 

91) we obtain 



n2 



1. Using relation and 

(l)^ ^lAi^-^i 



m 



T(A;)^(A;} 



'S'2'2(^l ^ ^) 



+ d{k) 



5*21 (fc, fci 



+ d{k 



(1)n ^l',2{kl \ k) 



B{k)\^ >= A{k)^{kl 



which is clearly satisfied if the coefficient of B{k)\Q > (unwanted term) vanishes, 
condition, which fixes ki^\ is given by 



Y[slS{kp,,k?) 



1, 



where we used the relation 



valid for the 5" matrix 



Aik) 



S'2'x (A;, A;') ^ Si'2{k',k) 
Si'2 {k, k') Si^2 (A^'i ^) 
). The eigenvalue is given by 
Sii{k^k\^) S2]2{ki\k) 



5']^'2(/c, k[ ^) Si^2{^i ^ A;) j=i 
provided (0)is satisfied. 

2. In this case the application of A{k) and D{k) in the ansatz (|48| ) give us 



+ 



:i-2/,(i) 



IlSliikp^,k), 



n2 



(1) Ji)^ 



Si I {k, k[ ^)Si i {k, k2 \ 



Si 2{k,ki ^)Si 2{k, k2 



(1) 



a(A;f^)S(A;)5(A;^'0l^ > 



and 



D{k)^{k[ 



(1) 



Si'2{k, ki ^)Si^2{^i \ ^2 ^) 
'S'2,'1 (A^) A;2 '*)<S'i'i {k2 \ k\ '^) 
Si'2{k, k2 ^) 31^2(^2 \ki "*) 

, A;) 5*2 '2(^:2 \ k) 



a{k^2^)B{k)B{k['')\Q > 



C2,2/,{1 
^2,2l'^l 



S*! 2 (M 
2(M 



(1 



5*1 '2 {ki 

3 1,2 i^2 



3 1,2 {^2 



, k)Si 2{k2 \ k) 
, A;) 5*2 '2(^:2 A;j; ■* 



'^2 y 



, A;)S'i'2(A;2 A;^ "* 
, A;) 5*2 '2 ( A: j A;2 



, A;)S'i'2(A;( A;2 "* 



d{k)^{kY 
d{k'i^)B{k)B{k^2^)\n > 

d{k^2^)B{k)B{k'i^)\n >, 
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where besides the relation and (^51) we have used (|43|) , which imply \E'(fcJ ,k 



(1) um^ 



2 



^'(/c2^'', /cj^''). From ( ^41) and (|55|) the condition (]49| ) gives us 



n 2 



cl.lfL n 2 C2,2.,(l) ,^ 

A(*) = n n + n «=) n leij^- 

j=i j=i '-'1,2 1"') "-j ) '=1 i=i '-'1,2 l^-j 5 "-J 



under the condition, which fixes k^^ and k2'\ 

ns;j(»=P.,4")= n i ngiiifc^.^i^') "=1.2- (57) 



In deriving this last expression the relation ( p2D was also used. 

General n^. The previous procedure can be iterated straightforwadly for arbitrary numbers 
n2 of second-class particles, which gives 



where 



(1) 



i+l) 



i=l 



— ^ A- ^(A;)^(A;, A;( \ . . . , fc^Ji, k\^i, 



i=l 



1 ^n2 ) ' 



(58) 



A(^)(A;) 
A(^)(A:) 



Ar(^) 



n2 



Si'i {k, k, 
i=i Si^2 {k, k\ ■*) 
7 02^2 1'^j ) '^1 



i=i «=i '-'1,2 l^j ) 

^1,2/,. ,„(1) 



(A)/, N _ 'S'2;i (A;, kj ) 



a(fcf)) n 



n2 



S'1'2 (A;, A;^ ) 
^1,2 {kj \ k) 



^1,2 {kj \ k) 



d{kf^] 



«=i(^j) .bi^2(% ,A;^ j 

n2 C2.2/7,(1) 7,(1)\ 

n '^2,2 y^l ) '^j J 
el, 2/7 ( 



(59) 



Then ^(A:f\...,A;«) is an eigenvector of T with eigenvalue A(/c) = A('^^(A;) + A*^'^^(A;), i. e. 



=1 2 \ / «=1 a=l '^12 ) ) 
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if the following conditions, which fix {ki^\ . . . , kn^} are satisfied 



A^\k)+A^\k)=0, z = l, 



,n2. 



(61) 



Using the relations (|52D, this last condition can be written as 

" ^a(fcp^,fcW) - si:M\kf^)sii{k^^\k(^y) 

IItXTT^ rTTTT= 11 



a = 1, ... ,712, 



(62) 



J=l ^1,1 1'^-Pj ? '^a J /3=l{7^a) ^1,2 I'^a 5 '^/3 J*-'2,2l'^/3 i J 

which concludes the diagonalization of T{k). 

Now let us return to our original problem of finding the eigenvalues of the Hamiltonian 
(^). The Bethe-ansatz equations will be obtained by inserting in (|38D the eigenvalues eval- 
uated at kj, i. e., A{kj), given in (|60D, with the condition (0). Taking into account that 
Si 2{k, k) = 0, we obtain 



-ikiN 



«2 



n(s,A':;(*.,*,))n||M''' 

1=1 a=l *Jl,2 ) ""O 



,n. 



(63) 



with the condition 



^ sl^{k,M^)) _ - s]'^{k^\k^^)s]:i{k^\kO:)) 

11 ..1.1/, , fl)N 11 



a= 1, . . . , ^2. 



(64) 



Writing explicitly the S-matrix elements ( p6D in (^) and ( [6^ ) we can state that the energies, 
in the sector with ni first-class particles and n2 second-class particles are given by 



where {kj,j 



E = -Y,{e.e'^^ + e+e-'^^ - 1), 
, n} are given by the solutions of 



^jfcj(Af-7l2(S2-Sl)) 



(_)n.-l]-[ 



112 

X n 



1=1 

e+(e 



i{ki+kj) 



^ikj 



e+ + e_e 



i{ki+kj) 



^iki 



(l)^ 



(i = l,•••,?^), 



(65) 



(66) 



Jkj{s2-si)n ( \n 



n 



\ni-l 



n2 g 

/3=i e+ + e_e^ " ^ " ' - e " 



(67) 
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It is interesting to observe that in the particular case where n2 = we obtain the Bethe- 
ansatz equations, recently derived [|1^], for the asymmetric diffusion problem with particles 
of size Si. Also the case si = S2 = 1 give us the corresponding Bethe ansatz equations for 
the standard problem of second class particles. The Bethe-ansatz equations for the fully 
asymmetric problem is obtained by setting e_|_ = 1 and e_ = 0. 



IV. CONCLUSIONS AND GENERALIZATIONS 

We obtained through the Bethe ansatz the exact solution of the problem of first- and 
second-class particles diffusing and interchanging positions on the lattice. We show that the 
solution can be derived in the general case where the particles have arbitrary sizes. 

Some extensions of our results can be made. A first generalization is the problem where 
M > 2 hierarchical ordered classes of particles (1st, 2nd, ...,Mth classes), with sizes Sj 
[i = 1,...,M), besides diffusing on the lattice, interchange positions in such way that 
particles of higher hierarchical classes see the lower ones in the same way as they see the 
holes (0 class). The allowed processes and the Hamiltonians are simple generalizations of 
relations (@)-(i) and @, respectively. The Hamiltonian obtained in a open lattice will be 
Uq{SU{M)) symmetric and in the case Si = S2 = • • • = -Sm = 1 if is related with the spin 
S = Y Sutherland model Jl^ . The model is certainly integrable and its eigenspectra will be 
given in terms of M-nested Bethe ansatz equations which are generalizations of (|66D-(|67|). 

A further quite interesting generalization of our model happens when we consider 
molecules in one or both classes with size s = 0. Molecules of size zero do not occupy 
space on the lattice, having no hard-core exclusion effect. Consequently we may have, at a 
given lattice point, an arbitrary number of them. The Bethe-ansatz solution presented in 
the previous section is extended directly in this case (the equations are the same) and the 
eigenenergies are given by fixing in (|65D-(|6^) the appropriate sizes of the molecules. It is 
interesting to remark that particles of second class with size S2 = 0, contrary to the case 
S2 > 1, where they "accelerate" the diffusion of the first class particles, now they "retard" 
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the diffusive motion of these particles. The quantum Hamiltonian in the cases where the 
particles have size zero are written in terms of spin s = oo quantum chains. 

The Bethe-ansatz equations in the case of the asymmetric diffusion, with particles of unit 
size [0,0, or with arbitrary size |T^, was used to obtain the finite-size corrections of the 



mass gap Gat of the associated quantum chain. The real part of these finite-size corrections 
are governed by the dynamical critical exponent z, i. e., 

Re(Gjv) ~ A^"'. (68) 
The calculation of the exponent z for the model presented in this paper, with particles of 



arbitrary sizes, is presently in progress |24 . 
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Figure Captions 



Figure 1 - Example of configurations of molecules with distinct sizes s in a lattice of size 

= 6. The coordinates of the molecules are denoted by the black squares. 
Figure 2 - Graphical representations of: a) the S matrix (|26|) and b) the Yang-Baxter 
equations (p^. 

Figure 3 - Graphical representation of the relation (0). The circles represents the particle's 
indices {Q} and {Q'}, and the lines the "wave numbers" {kp^}. The variables represented 
in the full circles are summed. 

Figure 4 - a) Graphical representation of the result of n iterations of relation (|3lD on the 
right of equation (^31). The diagramatic representation of Fig. 3 was used, b) A graphical 
representation of equation (|36D . 

Figure 5 - Graphical representation of the transfer matrix T{k), of the inhomogeneous 6- 
vertex model, in a periodic lattice of size n. 

Figure 6 - Graphical representation of the relation ( ^0] ) between the S matrix and the 
monodromy matrix Ai. 

Figure 7 - Graphical representation of the relations (^) . 
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